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Abstract 

Entanglement degradation caused by the Unruh effect is discussed for the tripartite GHZ or W states 
constructed by modes of a non-interacting quantum field viewed by one incrtial observer and two uniformly 
accelerated observers. For fermionic states, the Unruh effect even for infinite accelerations cannot completely 
remove the entanglement. However, for the bosonic states, the situation is different and the entanglement 
vanishes asymptotically. Also, the entanglement is studied for the bipartite subsystems. While for the GHZ 
states all the bipartite subsystems are identically disentangled, for the W states the bipartite subsystems 
are somewhat entangled, though, this entanglement can be removed for appropriately accelerated observers. 
Interestingly, logarithmic negativity as a measure for determining the entanglement of one part of the system 
relative to the other two parts, is not generally the same for different parts. This means that we encounter 
tripartite systems where each part is differently entangled to the other two parts. 
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I. INTRODUCTION 



The phenomenon of entanglement has a central importance in the quantum information science, 
and has emerged as a fundamental resource in quantum communication, quantum cryptography, 
quantum teleportation and quantum computation [1, 2]. Recently, much attention has been given 
to relativistic effects in the context of quantum information theory. Understanding entanglement 
in a relativistic setting is important both for providing a more complete framework for theoretical 
considerations and for practical situations such as the implementation of quantum computation 
tasks performed by observers in arbitrary relative motion. So, the relativistic quantum information 
theory may become a necessary theory in the near future. Relativistic quantum information in 
inertial frames has already been studied [3-8]. Peres et al demonstrated that the spin of an 
electron is not covariant under Lorentz transformation [3] . Alsing and Milburn [4] studied the effect 
of Lorentz transformation on maximally spin-entangled Bell states in momentum eigenstates and 
Gingrich and Adami [5] derived a general transformation rule for the spin-momentum entanglement 
of two qubits. Now, it is well known that for different observers in uniform relative motion the 
total amount of entanglement is the same in all inertial frames, although they don't agree on the 
amount of entanglement among various degree of freedom of the system. 

Also, more recently, quantum entanglement has been studied in relativistic non-inertial frames 
[9-17]. Alsing and Milburn extended the argument to a situation where one observer is uniformly 
accelerated [9]. Fuentes-Schuler and Mann calculated the entanglement between two free modes 
of a bosonic field, as seen by an inertial observer detecting one mode and a uniformly accelerated 
observer detecting the other mode [10]. Alsing et al did this calculation for two free modes of 
Dirac field [11]. Pan and Jing discussed the degradation of entanglement for two non-maximally 
entangled free modes of scalar and Dirac fields [12]. Martin-Martinez and Leon discussed the 
behavior of quantum and classical correlations among all the different spatial-temporal regions of 
a spacetime with an event horizon for both fermionic and bosonic fields [14]. In another work they 
relaxed the the single-mode approximation and analyzed bipartite entanglement degradation due 
to Unruh effect by introducing an arbitrary number of accessible modes [15]. Also Bruschi et al 
addressed the validity of the single-mode approximation that is commonly invoked in the analysis 
of quantum entanglement in noninertial frames. They showed that the single-mode approximation 
is not valid for arbitrary states, finding corrections to previous studies beyond such approximations 
in the bosonic and fermionic cases [16]. 

In this paper, we attend to tripartite systems which are constructed with modes of a non- 
interacting fermionic or bosonic field. For our purpose, we constraint our argument to the single- 
mode approximation [9-14]. It is shown that as a universality principle going beyond this ap- 
proximation does not modify the way in which Unruh decoherence affects fermionic and bosonic 
entanglements. Unruh effect is independent of the number of field modes considered in the prob- 
lem and statistics is the ruler of this process [16]. In tripartite discrete systems, two classes of 
genuine tripartite entanglement have been discovered, namely, the Greenberger-Horne-Zielinger 
(GHZ) class [22, 23] and the W class [24, 25]. These two different types of entanglement are not 
equivalent and cannot be converted to each other by local unitary operations combined with clas- 
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sical communication [19]. The entanglement in the W state is robust against the loss of one qubit, 
while the GHZ state is reduced to a product of two qubits. According to the geometric measure 
of entanglement, the W state has higher entanglement than the GHZ state does [26]. Methods are 
proposed for generation and observation of GHZ or W type entanglements [27, 28]. 

Here, we consider three observers, two non-inertial observers with different uniform accelerations 
and one inertial observer. They initially share a state built by modes of the fermionic or bosonic 
field, as viewed in the Minkowski coordinates. The state is chosen to be entangled as the GHZ 
state or the W state. Then switching to the Rindler coordinates for the accelerated observers, we 
investigate the change of entanglement caused by the Unruh effect [18]. We use the logarithmic 
negativity as a measure to evaluate the entanglement. This measure can be used to evaluate 
the entanglement of one part relative to the other parts of a system with arbitrary dimension 
[20] . In a similar work Hwang et al have discussed tripartite entanglement in a non-inertial frame 
[17]. However, in their problem only one of the three observers is accelerating and only bosonic 
entanglements are discussed. 

This paper is organized as follows. In Sec. 2 we review the basic formalism in brief and introduce 
the Bogoliubov transformation both for fermionic and bosonic states. In Sec. 3 we discuss the 
Unruh effect for the fermionic GHZ or W entanglements. We use logarithmic negativity as a 
measure for the entanglements. The effect of Unruh temperature on residual entanglement of 
bipartite entanglements is also discussed. In Sec. 4 similar to Sec. 3 we address the entanglement 
degradation for Bosonic GHZ or W entanglements. We demonstrate our results by appropriate 2 
or 3 dimensional graphs. Finally in Sec. 5 we present the conclusions. 

II. BASIC FORMALISM 

Prom an inertial perspective, Minkowski coordinates are the most suitable coordinates for dis- 
cussing the problem. However, for a uniformly accelerated observer the Rindler coordinates are 
appropriate. Two different sets of Rindler coordinates are required for covering the Minkowski 
space. These two sets, define two causally disconnected regions. Therefore, the accelerated ob- 
server has access only to one of the regions and so he must trace over the states in the inaccessible 
region. This process leads to an information loss and is known as the Unruh effect. During the 
Unruh effect the pure state of one region evolves to a mixed state and so the entanglement de- 
creases. This effect can be considered as a local non-unitarian transformation, since acceleration 
affect every mod locally. The Unruh effect implies that a vacuum state in an inertial frame is seen 
as a thermal state by an accelerated observer, where the corresponding temperature depends on 
the value of acceleration. 

Consider a uniformly accelerated observer moving with a proper acceleration a in Minkowski 
(t, z) plane. Two different sets of Rindler coordinates (r, ^) defined as 

at = e°'^ sinh or az = e"^ cosh or region I (1) 

at = — e"^ sinh ar az = — e"^ cos ar region II 



3 



are necessary for covering Minkowski spacetime [21]. These coordinate transformations define 
two Rindler regions I and II that are causally disconnected. In Minkowski coordinates, let the 
operators {Ak-,A\.) be annihilation and creation operators for the positive energy solutions and 
{Bk,Bl.) be annihilation and creation operators for the negative energy solutions. Here, A; is a 
notational shorthand for the wave vector k, which labels the modes. Also, in region I, let us denote 
{C\,C^i}) as the annihilation and creation operators for particles and as the annihilation 

and creation operators for antiparticles. The corresponding particle and antiparticle operators in 
region II are denoted as (C^^C^^^) and {'D^,V^^^). 

As we know the relationship between the Minkowski and Rindler creation and annihilation 
operators is given by the Bogoliubov transformation. Thus, the Bogoliubov transformation for a 
fermionic field is 



(2) 



where tan it = e '^^'^/'^ with lo as the particle frequency and c as the speed of light. Also, for a 
bosonic field we have 
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where tanhr = g,-'^^^/'^. Note that < tt < 7r/4 and < r < oo, as the proper acceleration a takes 
its full range, i.e. < a < oo, 

Using the above Bogoliubov transformations we can express the vacuum Minkowski state |Ofc)^ 
and the first excited state in terms of Rindler states. For the fermionic field we can write 



\^k)X[ = cosn|OA;)j'"|0_fc)ji + sinn|lA;)i'"|l-fc)n 



\^k)\i - |lfc)l'|0-A;)n- 



(4) 



Also, for the bosonic field we have 

^ oo 



n=0 

oo 
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\^ V tanh" rVn+ 1| (n + l)fe)j+ 



cosh r 



n-k)^l- 



(5) 



n=0 



In these transformations \n\^Y ^'^d indicate the Rindler-region-I particle state and Rindler- 

region-II antiparticle state, respectively. 
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FIG. 1: Rindlcr spacetime diagram for an inertial observer Alice and two accelerated observers Rob and 
Steven. The thick black line shows the word line of Alice (A). A uniformly accelerated observer Rob (R) 
with acceleration ai travels on a hyperbolic world line constrained to region I, and a fictitious observer 
anti-Rob (R) travels on the corresponding hyperbola in region II given by the negative of Rob's coordinates 
(two thick red hyperbola). The horizons of R (R) are lines of r = ±00 which Alice (A) will cross them 
at finite Minkowski times corresponding to T-ij^ (two dashed blue lines). Another uniformly accelerated 
observer Steven (S) travels on another hyperbola with an acceleration 02 constrained to region I', and a 
fictitious observer anti-Steven (S) travels on the corresponding hyperbola in region II' given by the negative 
of Steven's coordinates (two thick green hyperbola). The horizons of S (S), are lines of r' = ±00 which Alice 
(A) will cross them at finite Minkowski times, are (two dashed cyan lines). 

III. FERMIONIC ENTANGLEMENTS 
A. The GHZ state 

Consider three observers, Alice, Rob and Steven, such that Alice is at rest in an inertial frame, 
but Rob and Steven are moving with uniform accelerations with respect to Alice. We assign 
accelerations ai and 02 to Steven and Rob, respectively. Fig. 1 shows the corresponding spacetime 
diagram. We see that at some points Alice's signals will no longer reach to Rob and Steven, 
however, Rob's and Steven's signals will always reach to Alice. Each observer carries a detector 
sensitive only to a single mode of a fermionic field, for Alice, /cr for Rob and for Steven. We 
suppose a GHZ entanglement for this tripartite system, as viewed in an inertial frame. Therefore, 
we use the Minkowski modes to construct a GHZ state for the system as 

|GHZ)^Rg = ^ (|0fcj + |0fc,)+|0fc,)+ + + . (6) 

where we have used this fact that in a fermionic field there are only two allowed states for each 
Minkowski mode. 

Now, we need to express the states jn^j^)^ and |nfcg)~^, in terms of Rindler states corresponding 
to Rob and Steven. Then using (4), we obtain 

|GHZ)ars = [|OfeA)'^(cosui|Ofcj,)j+|0„fcp)j"^ -Fsinni|lfcj,)J^|l_fej,)j"^)(cosu2|0fcg)^|0_fe5)j"j, 

+ sinn2|U3)+|l-fcs)rr) + |UA) + (|Us)i+|0-fcs)„)(|lfcs)it|0_fcs)jj,)] , (7) 
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where \nkj^)f (or \nkg)y) is the Rindler-region-I-particle mode for Rob (or Steven), |n_fcj^)jj (or 
|n_fcg)jj,) is Rindler-region-II-antiparticle mode for Rob (or Steven), and tanuj = exp(— vrtjc/aj) 
with i G {1,2}. 

Using the basis {\nkj+\nk^}+\nks)+} that is {|000), |001), |010), |100), |011), |101), |110), |111)}, 
we can obtain the density matrix for the GHZ state (6) as 

Pars = ^ (Si (g) Si (g) Si - Si (g) S2 O S2 - S2 Si S2 - S2 O S2 ® Si 

o 

+S3 (g) S3 (g) / + S3 (g) / (g S3 + / (g) S3 (g) S3 + / (g) / (g) /) (8) 

where / is the unit 2x2 matrix and {Si, S2, S3} are the Pauh matrices. Notice that here and in 
the fohowing, the matrices in each tensor product term are placed according to Ahce-Rob-Steven 
order. The pure state (8) describes a tripartite system. On the other hand, the density matrix for 
the state (9) is obtained as Pa,i,i',ii,ii' = |GHZ)ars(GHZ|, which is pure and describes a five-partite 
system. However, as Fig. 1 represents, whole of spacetime is accessible only for the inertial observer 
Alice, and the accelerating observer Rob (or Steven) has only access to one region say I (or I'). So, 
we must trace over the states belong to regions H and H'. Doing so, we get the following tripartite 
density matrix 

Pa 1 1' = - [costii costi2(Si (g) Si (g) Si - Si (g) S2 (gi S2 - S2 g) Si g) S2 
' ' 8 

-S2 S2 Si) + ^(cos2ni cos2u2 - 1)(S3 S3 ® S3) + cos^ mi(S3 (g S3 /) 

+ cos^ U2(S3 0/0 S3) + ^(cos2ni cos 2^2 + 1)(/^ S3 ® S3) - sin^tii(/® S3 ® /) 
-sin2'U2(/(g/(g)S3) + /(g/(g)/] (9) 

where we have employed the basis {\'nkj^)'^\nkji)f\nkg)^,}- While the state (8) is pure and maximally 
entangled, the state (9) is not pure and, as we will verify, its entanglement is degraded. This 
entanglement degradation is essentially justified by the Unruh effect. 

1. A-RS, R-AS, S-AR entanglements 

To quantify the entanglement in a tripartite system, different measures have been introduced. 
Here, regarding the dimension of the density matrix (9) and that the state is not pure, we consider 
the logarithmic negativity which describes the entanglement of one part of the system relative to 
the other parts. For instance, the logarithmic negativity of Alice part relative to the other two 
parts is defined as A/a-rs = Z]jlog2l'^i| where Aj denotes the eigenvalues of /o^jj/ which is the 
partial transposition of 1 1' with respect to Alice. Similarly, we can calculate A/r-as and As-ar 
by finding the eigenvalues of p^jj/ (partially transposed density matrix with respect to Rob) 
and J ~, (partially transposed density matrix with respect to Steven) . Logarithmic Negativity 
vanishes unless some negative eigenvalues are present. Let N denote the negative eigenvalue, then 
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we can write the logarithmic negativity also as 



AA = log2(l-2Af) (10) 

We can readily obtain the required partially transposed matrices from (9) by noting that after 
a transposition, Si and S3 do not change but S2 — > — S2. We have 

/Ot T J, = £1 - cosiii cosu2(Si (g) S2 <8) S2 - S2 (8) Si (g) S2 - S2 <8) S2 <8) Si), (11) 



and 



Paii' = Q - -^cosui cosn2(-Si ® S2 ® S2 + S2 ® Si ® S2 - S2 S2 Si), (12) 



J p = ^ — - cosni cosn2(— Si (g) S2 (g S2 — S2 (g Si (g) S2 + S2 (g S2 (g Si), (13) 



where 



cos til costi2(Si (g Si (g Si) + -(cos2ni cos 2^2 — 1)(S3 (g S3 (g S3) 



(14) 



+ cos^ til (S3 (g S3 (g /) + cos^ ^2(^3 (g / (g S3) + -^{cos 2ui cos 2U2 + 1)(/ (g S3 (g S3) 

,2 „, /T ^ T\ „;„2 



sin^ ui{I S3 /) - sin^ U2{I / (g S3) + / ® / (g / . 



The eigenvalues of each of these partially transposed matrices can be obtained explicitly, however, 
we need only negative eigenvalues. It turns out that the negative eigenvalues for (11), (12) and 
(13) are 



-^A-RS = ^ sin^ ui sin^ U2 — - V sin^ ui sin^ U2 + 4 cos^ ui cos^ U2, 



(15) 



-^R-AS = ^ sin^ ui cos^ U2 — - \/ sin^ ui cos^ ?i2 + 4 cos^ ui cos^ U2 , (16) 

and 

-^S-AR = ^ sin^ U2 cos^ ui — - \l sin^ U2 cos^ Mi + 4 cos^ U2 cos^ u\ , (17) 

respectively. For investigating the entanglement, we substitute these eigenvalues in (10) to obtain 
the corresponding logarithmic negativity as functions of u\ and U2- Recall that in the present 
case < < ^ corresponding to < < 00 where Oj indicates the proper acceleration of Rob 
or Steven. We have plotted all three surfaces A/r,_aS) A/s-ar and A/a-rs in Fig. 2(a). Each 
surface shows a uniform decreasing of logarithmic negativity as accelerations increase. It must be 
noted that even for infinite accelerations, that is, ui,U2 — ?• f , each logarithmic negativity does not 
vanish. In other words, like bipartite entanglement, the Unruh effect doesn't completely destroy 
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FIG. 2: Logarithmic negativity for the fermionic GHZ state versus accelerations ui and U2- (a) The two 
lower surfaces that intersect at ui = U2, represent A/r-as and A/s-ar- The upper surface corresponds to 
A/a-rs- (b) A section of figure (a) for a given U2 = 37r/16. Even for infinite accelerations or ui, U2 — >■ j, each 
logarithmic negativity does not vanish, i.e., the Unruh effect does not completely destroy the entanglement 
in this case. 

the entanglement. We see that the surface of A/a-rs covers the other two surfaces. To illustrate 
the situation a section of Fig. 2(a) for a given U2 = Sn/IG is plotted in Fig. 2(b). It is remarkable 
that these entanglements are not generally equal. This means that in the considered tripartite 
system each part is differently entangled to the other two parts; the inertial Alice part is mostly 
entangled. Of course, the surfaces of A/r-as and A/s-ar intersect at ui = U2, i.e., when Rob and 
Steven have the same acceleration, have the same entanglement, expectedly. 



2. Entanglement of bipartite subsystems 

In this subsection, we attend to the entanglement of bipartite subsystems AR, AS and RS. 
To do this, we must trace out one of the parts of the tripartite system ARS described by the 
density matrix (9). We can readily trace over Alice, Rob or Steven states by noting that the Pauli 
matrices are traceless. For instance, upon tracing over the Alice states, only the tensor product 
terms remain where the first matrix is the unity matrix I. Then, we obtain the following reduced 
density matrix for RS subsystem 



1 



^ (cos 2^1 cos 2^2 + 1)(S3 (g) Ss) -sin^ ui{T,3 I) - sin^ n2(/ (8) S3) + / (g) / 



(18) 



which is diagonal. In the same manner we see that the reduced density matrices for AR and AS 
subsystems are also diagonal. Thus, each of these bipartite subsystems is disentangled. This means 
that the GHZ character of the state is preserved under the Unruh effect, that is, tracing over any 
part of a GHZ state, leads to a disentangled bipartite subsystem [19]. 
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B. The W state 



In this subsection we assume a W entanglement for the tripartite system ARS as viewed in an 
inertial frame. Therefor, we use the Minkowski modes to construct a W state as 



|W)fRS = ^ (|1..>+|0.,)+|0,,)+ + |0,,)+|1,,)+|0,,)+ + |0,J+|0,,)+|1,,)+) . (19) 

The corresponding density matrix for this state is obtained as 

Pars = ^ [2(Si ^ Si ® S3 + S2 ® S2 ® S3 + Si ® S3 ® Si + S2 ® S3 ® S2 (20) 
+S3 ® Si (g) Si + S3 (g) S2 S2 + Si / Si + Si Si / + S2 S2 / 
+S2 ® / S2 + / ® Si ® Si + / S2 S2) - 3(S3 O S3 S3) - S3 ® / S3 
-/ S3 S3 - S3 «) S3 / + S3 ® / / + / ® S3 ® / + / ® / S3 + 3(7 / ® /)] . 

Then we apply the Bogoliubov transformation (4) for the accelerating observers Rob and Steven 
and rewrite the state (19) in Rindler coordinates which leads to a density matrix describing a 
five-partite system. After tracing over the causally disconnected regions II and II', we reach to 

Pa,i,i' = ^ [2 cos Ml cos2ti2(Si (g) Si (g) S3 + S2 S2 S3) + 2cosn2(Si (g)/(g)Si + S2®/(8) S2) 
+2cosni(Si (g Si (g / + S2 (g S2 (g /) + 2 cos 2ni cosn2(Si (g) S3 (g) Si + S2 (8) S3 (g) S2) 
+2 cos ui cos ti2(S3 (g Si (g Si + S3 (g S2 (g S2 + / Si Si + / (g S2 S2) 

— (cos 2ni cos 21^2 ~1~ COS 2l/i -\- cos 2li2)(S3 (g S3 {g S3) — S3 {g S3 {g / — S3 (g / <g S3 -j- S3 {g / g) 7" 

+(cos 2ni cos 2u2 — cos 2ni — cos2n2)(/ (g S3 (g S3) + (2 cos 2tii — 1)(/ (g S3 g) /) 

+(2 cos 2n2 - 1)(/ g) / <g S3) + 3(1 g) / g) /)] (21) 



1. A-RS, R-AS and S-AR entanglements 

Again we want to use the logarithmic negativity given by (10) for evaluating the entanglement of 
the state (21). To do this we first obtain the partially transposed matrices PaiF'^aTi' ^'^'^ '^ai?' 
Then we must calculate the negative eigenvalues for them and substitute in (10). However, these 
calculations are lengthy and here we only present the results by plotting the obtained logarithmic 
negativity in Fig. 3. These surfaces show the logarithmic negativity in term of accelerations ui and 
U2- We see that the surface of A/a-rs lies below of the surfaces of A/r,_as and A/s-ar, and is covered 
by them when ui and U2 take their full range. Fig. (3(c)) represent the situation more clearly. As 
is seen, again in the considered system each part is entangled differently. Expectedly, the surfaces 
A/r-as and A/s-AR intersect at ui = U2. However, it is interesting to note that the surface A/a-rs 
intersects the A/r-as and As-ar at definite values of ui and U2-, i.e., in the considered W state, for 
definite accelerations the entanglement of part R or S can be equal to the entanglement of inertial 
part A. Note that similar to the GHZ case discussed in the previous subsection, some degree of 
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(a) 



(b) 




(c) 

FIG. 3: Logarithmic negativity for the fermionic W state versus accelerations ui and U2. (a) The surface 
A/a-rs in terms of ui and U2- (b) A/a^rs, A/r-as and A/s-ar surfaces all together. Again Ar_as and 
A/s-AR surfaces intersect along ui = U2- The surface of A/a-rs is not seen here , because, this lies below 
the other two surfaces, (c) A section of the figure (b) for a given U2 = tt/S. The surface A/a-rs intersects 
Ar,-as or A/s-AR at determined values of ui and U2- Like the GHZ case, even for infinite values of oi and 
a2, the logarithmic negativity generally is nonzero. 

W entanglement is preserved as ai and 02 go to infinity. This result is obtained also for fermionic 
bipartite entanglement and has been proven to be universal for fermionic fields [15]. 

2. Entanglement of bipartite subsystems 

Let us study the entanglement of bipartite subsystems AR, AS and RS, by tracing over the 
states of one of the parts of the tripartite system described by (21). We denote the resulting 
density matrices as /3jj/, pA,i and />ai') when we trace over Alice, Steven and Rob, respectively. 
We readily obtain from (21) 

Pi I' = [2 ui cos ti2(Si ig) Si + S2 (X" S2) + (cos 2ui cos 2u2 — cos 2ui — cos 2n2)(S3 S3) 

+(2cos2ni - 1)(S3 0/) + (2cos2n2 - 1)(/ S3) + 3(1 ® /)] , (22) 
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PA,i = ^[2cosni(Si(8)Si + S2i^S2)-S3(g)S3 + S3(g)/ (23) 
+ (2 cos lux - ® S3) + 3(7 ® I)\ , 



and 



Pa,!' = 3^ [2 cos ti2(Si ® Si + S2 1:2) - S3 ® S3 + S3 / (24) 
+ (2 cos 2n2 - 1) (/ S3) + 3(1 /)] , 

In contrast to the GHZ case, these bipartite subsystems are not disentangled. We can still use the 
logarithmic negativity (10) for calculating these bipartite entanglements. To do this we must first 
find the negative eigenvalues for partially transposed matrices Pjj/, and corresponding to 
the above matrices. We obtain 

A'rs = (cos^ ui + cos^ 7x2) H — (cos ui cos ^2)^ 

2 3 3 

— Y^9 — 12(cos2 Ml + cos^ U2) + 12(cos u\ cos u-i)^ + 4(cos^ u\ + cos^ ^2), (25) 
6 



1 1 



iVAR = 7: - - Vl + 4cos4n2, (26) 
b b 



1 1 



Nk^ = --- Vl + 4cos4ui, (27) 
b 

which as substituted in (10), give the corresponding logarithmic negativity. It is remarkable that 
the negativity iY^s in (25) and consequently the logarithmic negativity A/rs vanish if u\ and U2 
satisfy the equation 



\/2sin?/i 
3" 



COSM2 = ^ 2 • (^^) 

V 2 — cos^ u\ 



Fig. 4(a) shows the behavior of A/rs in terms of u\ and U2- We see that if the accelerations of Rob 
and Steven satisfy (28), the entanglement between them will be removed, completely. This is a 
remarkable result and it seems that this contradicts the general behavior of fermionic entanglements 
under the Unruh effect. However, A/rs represents the residual entanglement of Rob and Steven 
parties after tracing over the Alice states. Thus, the entanglement level for this system is lower 
than the entanglement of the whole tripartite system, as Fig. 4(a) shows. Moreover, both Rob 
and Steven are accelerated observers and so the rate of entanglement degradation is such that the 
entanglement descends to zero for finite accelerations. 

In Fig. (4(b)) we have plotted A/ar (A/as) which depends only on u\ {u^). Since the entangle- 
ment is destroyed upon tracing over Steven (Rob), the curve starts with 0.5 at ui = {u2 = 0), as 
for RS subsystem. But, here only one of the observers accelerates and the entanglement degradation 
is not enough for vanishing the entanglement at a finite or even infinite acceleration. 
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FIG. 4: The logarithmic negativity versus ui and U2 for bipartite subsystems when the tripartite system 
is in the fermionic W state, (a) A/rs surface which represents the entanglement of Rob and Steven states 
in their first regions I and l'. As is seen, A/rs equals 0.5 at ui = U2 = 0, then reaches to zero at a curve 
defined by (28). (b) A/ar (A/as) depends only on ui (1x2). The curve starts with 0.5 at ui = {u2 = 0), 
then descends to a non-zero asymptotic value. 

It must be noted that after tracing out any part of the tripartite state (21) we obtain a bipartite 
subsystem having some residual entanglement, in an apparent contrast to the GHZ case. One 
may say that the Unruh effect does not change the class of the W state (19). Notice that a W 
state remains entangled after tracing out one of its parts. However, it is remarkable that for RS 
subsystem the residual entanglement can completely be removed for appropriate accelerations of 
Rob and Steven. 

IV. BOSONIC ENTANGLEMENTS 

In this section we are going to discuss the effect of Unruh temperature on the bosonic tripar- 
tite GHZ and W entanglements. Because of the nature of bosonic states and hence the form of 
corresponding Bogoliubov transformation (5), our calculations will be more complicated than the 
calculations for the GHZ case discussed in the before section. We begin with the bosonic GHZ 
state. 

A. The GHZ state 

Again we consider the GHZ state (6) which is written in terms of Minkowski modes. In the 
present argument we expand Minkowski states |Ofcj^)^ and |Ofcg)^ in terms of Rindler bosonic states 
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for Rob and Steven, using the Bogoliubov transformation (5). Then we can write 

|GHZ)ars = ^ |0)(^^^f;tanh"ri|n)i|n)„^ |^_L_ tanh^ r2|m)Hn)n' j (29) 

+ |1) f y^tanh" riVri + l|n+ l)i|n)ii ) 

ycosh n ) 

X ( \ — V tanh"" r2\/mTT|m + |m)n' ) • (30) 

ycosh r2^^ ) 

The corresponding density operator p = |GHZ)ars(GHZ| contains five partitions, however, as 
before we must trace out the causally disconnected regions II and II'. Then we reach to an infinite 
dimensional density matrix 

Pa 1 1' = 7T~ n — tanh^" ri tanh^™ r2 Pnm , (31) 

2 cosh ri cosh ro — 

^ ^ n,m 



where 

in \/n I , VnTT V^f^^^. , ^/n+l ^/m+l 

Pnm = \0,n,m){0,n,m\ -\ — 10, n, m)(l, n + 1, m + 1| H — 

cosh ri cosh r2 cosh ri cosh r2 

,ra+lm+l 

X |1, n + 1, m + 1)(0, n, m| H 5 2 — I-*-' n + 1,171 + l)\l,n + l,m + 1). 

cosh ri cosh r2 

1. A-RS, R-AS, S-AR entanglements 

In order to quantify the entanglement of the ARS system described by (31), we invoke the 
logarithmic negativity introduced in (10). First we must calculate the partially transposed density 
matrices Pah'' '''aTi' ^^"^ PaiV' "^^^^^ matrices have infinite dimensions, however they are block 
diagonal matrices. So, we encounter only square matrices located at each block. For instance, the 
(n, m) block of the Alice partially transposed density matrix, is obtained as 



X / tanh^ n tanh2 ra v^v^\ 

\ _ ^ ^ coshri coshr2 1 (•V)\ 

\ coshri coshr2 sinh^ ri sinh^ r2 / 

There are similar expressions for partially transposed density matrices for Rob and Steven. Let 
N^m be the negative eigenvalue of each block, then, the negative eigenvalue for the whole matrix 
is = '^^rn=Q-^rim which Can be used in (10) to get the logarithmic negativity. We obtain the 
negative eigenvalue for the block (32) as 



A^, 



A pq tanh^*^ ri tanh^"V2 / ,9 ,9 nm 
= 5 ^ tanh^ n tanh^ r2 + 



4 cosh^ ri cosh^ r2 \ sinh^ ri sinh^ r2 



2 

sinh'^ ri sinh'^ r2 / cosh'^ ri cosh"^ r2 



tanh^ ri tanh^ r2 + ^ ^ — 2 — ) ^ 4(n + m + 1) | ^ ^^^^ 
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(a) (b) 

FIG. 5: (a) A/a-rs, A/r-as and A/s-ar for the bosonic GHZ state. The surfaces A/r,-as and A/s-ar have an 
intersection at ri = r2 and the surface A/a-rs covers them. Here the entanglement asymptoticaUy vanishes, 
(b) A section of figure (a) for a given r2 = 1. 

Also, for the S-AR system, the negative eigenvalue for the (n, m) block is obtained as: 

AR tanh^"ritanh^''"r2 / ,9 (n + l)m 

4 cosn ri cosh r2 \ smh ri smh r2 



- / f tanh^ r + (^ + 1)"^ A ^ ^ 4(n + 1) 
y V sinh^ ri sinh^ r2 / cosh^ ri cosh^ r2 



It turns out that the negative eigenvalue for the (n, m) block of the R-AS system, can be obtained 
by interchanging r2 and ri in Eq.(34). 

Now, the logarithmic negativity A/a-rs becomes 

Ma-rs = log2 ( 1 - 2 f; N^™^' j , (35) 

\ n,m=0 j 

with similar expressions for A/s-ar and A/r_as- These functions are plotted versus r\ and r2 in 
Fig. 5(a). As the figure shows, the surfaces A/r-as and A/s-ar have an intersection at r\ = r2 
and the surface A/a-rs covers them, as for the fermionic GHZ case. A section of these surfaces for 
a given r2 = 1 is plotted in Fig. 5(b). The distinction between the curves again implies that in 
the considered tripartite system each part is differently entangled to the other parts. In contrast 
to the fermionic GHZ case, the logarithmic negativity in the present case asymptotically vanishes. 

2. Entanglement of bipartite subsystems 

Entanglement of bipartite subsystems AR, AS and RS can be obtained for each case by taking 
trace over the states in the otherwise subsystem. After doing some manipulations it turns out 
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that all the resulting bipartite density matrices are diagonal and so, we conclude that there is no 
entangled bipartite subsystem. This resembles the fermionic GHZ state. 



B. THE W state 

Now, let us apply the Bogoliubov transformation (5) in the W state (19) written in terms of 
Minkowski modes. Then we reach to the following state. 



|W)ars = ^ 



+ |0)a/ I Vtanh"riVn+ l|n + l)i|n)ii I f \ Vtanh'^r2|m)i/|m)n' ) (36) 

\^cosh2ri^^^ y\^coshr2^ J 

+ |0)a/ I T y^tanh"ri|n)i|n)n 1 ( 4 — "S^ tanh"" r2\/mTT|m + l)i'|m)ii/ 

\ cosh ri ^ / \ cosh^ r2 "H, 

\ n=0 / \ m=0 / 

Again tracing out the causally disconnected regions II and II', we reach to the following density 
matrix 

Pah' = z 72 72 — tanh^'' n tanh^™ r2 Pnm (37) 

dcosh ri cosh r2 — ^ 

n,m=0 

where 



\lnm){lnm\ + ~^ -|lnm)(0, n + 1, ml + H.C^ (38) 

\ cosh ri / 

+ ( '^"^^'^^ \lnm){0,n,m + 1| + H.C^ + ^ \^ |0,n + l,m)(0,n + l,m\ 
\ cosh r2 / cosh ri 



+ ( ^'^^^^"^^^ \0,n + l,m){0,n,m + l\ + H.c)] + |0,n,?n+l) (0,n,m+ 1|. 

cosh ri cosh r2 / cosh r2 



By inspection, we realize that the required partially transposed density matrices deduced from 
(37), are not block diagonal. So the calculation of logarithmic negativity for these density matrices 
by the trick of the previous subsection is impossible and we encounter a complicated problem that 
can be tackled by a numerical procedure. However, we do not follow this here and content ourselves 
with an approximation valid only for small ri and r2- Thus, we can consider the summation (37) 
up to m = n = 1, which leads to an 18 x 18 matrix. Then, we see that the behavior of the 
entanglements in the present case for small accelerations, is similar to what is shown in Fig. 5 for 
the GHZ state. 
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1. Entanglement of bipartite subsystems 



Let us trace out the Alice part of the state (37). Then, we reach to the following density matrix 
for the RS subsystem 

^ oo 

PiX = n 72 71— Yl t^nh^" ^1 tanh^™ /3„„^ (39) 

3 cosh ri cosh ro — 

n,m=0 

where 

Pnm = \nm){nm\-\ — n + 1, m)(n + 1, ?n H l,n + l)(l,m + l 

cosh ri cosh ri cosh r2 

m + 1 

H 2 — m + l)(n, m + 1|. (40) 

cosh r2 

After taking the partial transpose on Rob states we reach to a block diagonal matrix with the 
following negative eigenvalue for the (m, n) block 



tRS 



tanh2"ritanh2'"r2 / ^ /, ,,2 , , (n + l){m + l) 

-2 a + 6- W(a + b)^-4a6 + 4 ^ 2 

D cosh ri cosh r2 \ U cosh ri cosh r2 



where 



n m 

a = 1 H 2 1 2 ' 

tanh ri cosh ri tanh r2 cosh r2 

,9 9 (n + 1) tanh^ r2 im + 1) tanh^ ri 

h = 2 tanh^ n tanh^ r2 + o + - — ^• 

cosh ri cosh r2 

The corresponding logarithmic negativity A/rs is plotted in Fig. 6(a). As the figure shows, RS 
subsystem becomes disentangled at certain finite values of ri and r2. To determine these values 
one should find the root of (41). 

Now, tracing over the Steven states, we reach to the following reduced density matrix for the 
AR subsystem 

^ 00 

= T2— Yl tanh^" npn (42) 

3 cosh ri 

n=0 

where 



Pn = \l,n){l,n\ + (^^^^\ln){0,n + l\+ H.C. 
\ cosh ri 

n + 1 

+ 9— |0,n + l)(0,n + l| + |0,n)(0,n|. (43) 

cosh ri 

After taking the partial transpose on the Alice states, we reach to a block diagonal matrix that its 
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(a) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 O.S 0.9 1 



(b) 



FIG. 6: The logarithmic negativity versus ri and r2 for bipartite subsystems when the tripartite system is in 
the bosonic W state. (a) The logarithmic negativity A/r,s vanishes for certain values of ri and r2 determined 
by the root of (41). (b) The logarithmic Negativity A/ar (A/as) vanishes at ri{r2) = ln(l + \/2) corresponding 
to a certain value of Rob (Steven) acceleration. 



nth block has the negativity 

^2n 



A'™ = |^(l + -^+tanh^n (44) 

Dcosh ri V smh ri 



1 + h tanh^ ri ) - 4 tanh^ ri + ^ 



sinh^ ri J cosh^ ri 

which depends only on ri. Then, the logarithmic negativity is obtained by A/ar = 
log2 (1-2 J2'^=i ^n^) ■ It must be noted that (44) has a zero at ri = ln(l + \/2), independent of n. 
Then we have negativity only for < ri < ln(l + \/2). The logarithmic negativity A/ar is plotted 
in Fig. 6(b). As the figure shows, the entanglement vanishes at a finite acceleration corresponding 
to ri = ln(l + \/2). It can be shown that, the reduced density matrix for the AS subsystem has 
the same form of (42), but ri is replaced by r2. So, the results indicated in Fig. 6(b) can also be 
considered for the AS subsystem. 

Again note that the entanglement level for each bipartite subsystem is lower than that of 
whole system and here the Unruh effect can completely remove the residual entanglements at 
finite accelerations. Thus, for the bosonic W entanglement, all the bipartite subsystems become 
disentangled at finite accelerations, in contrast to the fermionic W states that its AR and AS 
subsystems never become disentangled (see Fig. 4). 



V. CONCLUSIONS 



Previously, in Ref. [15], the bosonic bipartite entanglement, and in Ref.[16], the fermionic 
bipartite entanglement, was discussed. However, in these works two observers were considered. 
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one inertial observer and one accelerated observer. However, our setting is more general and so we 
obtain some new results that are special to tripartite systems and distinguish this work from the 
previous works. 

In this work we considered the degradation of entanglement in tripartite entangled states caused 
by the Unruh effect. In particular, we considered two significant classes of tripartite systems namely 
GHZ and W states. These entangled states were built by three free modes of bosonic or fermionic 
quantum fields. One of these modes was observed by an inertial observer Alice and the other 
two modes were observed by uniformly accelerated observers Rob and Steven. This leads to the 
detection of thermal radiation by accelerating observers, which generally degrades the entanglement 
in the system. We showed that the Unruh effect, even for infinite accelerations, cannot completely 
remove the entanglement in the fermionic GHZ and W states. On the other hand for the bosonic 
states, we showed that the entanglement rapidly drops and is erased for large values of accelerations. 

We used the logarithmic negativity as a measure for these tripartite entanglements, and inter- 
estingly, the logarithmic negativity was not generally the same for different parts of the system. 
This means that we encounter tripartite systems where each part is differently entangled to the 
other two parts. For instance in the fermionic or bosonic GHZ state, the Alice part is mostly 
entangled to the Rob and Steven parts, for all accelerations. But for W states this depends on the 
accelerations. Of course, for determined accelerations it is possible that the entanglement be the 
same for two parts of the system. 

We also discussed the degradation of entanglement for bipartite subsystems. Both for fermionic 
and bosonic GHZ states, tracing over each part of the system leaves a disentangled bipartite 
subsystem. However, tracing out any part of the fermionic or bosonic W state leads to a bipartite 
system having some accelerated-dependent entanglement. It was deduced that for the fermionic 
W state, if the Alice part is traced out, the remaining entanglement can vanish for certain finite 
accelerations. But, if Rob or Steven part is traced out, the remaining entanglement will decrease 
to nonzero values, asymptotically. For the bosonic W state, we showed that all the bipartite 
entanglements can vanish for determined finite accelerations. 
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